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Abstract 



Oh; 

^ I We systematically derive a linear quantum collision operator for the spinorial Wigner transport 

r^ ' equation from the dynamics of a composite quantum system. For suitable two particle interaction 
potentials, the particular matrix form of the collision operator describes spin decoherence or even 

O |. spin depolarization as well as relaxation towards a certain momentum distribution in the long time 

rj ' limit. It is demonstrated that in the semiclassical limit the spinorial Wigner equation gives rise to 

^ ' several semiclassical spin-transport models. As an example, we derive the Bloch equations as well 
as the spinorial Boltzmann equation, which in turn gives rise to spin drift-diffusion models which 

^ I are increasingly used to describe spin-polarized transport in spintronic devices. The presented 

t^^ ' derivation allows to systematically incorporate Born-Markov as well as quantum corrections into 

■^ ■ these models. 
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I. INTRODUCTION 

The modeling of transport phenomena in electronic devices is one of the major challenges 
in modern solid state physics. While in most physical applications a full quantum mechanical 
treatment by means of the Schrodinger equation or the von Neumann equation is far too 
complex, it is a beneficial and legitimate approach to employ effective models. Well-known 
and prosperous examples are the drift-diffusion equations to treat systems in local thermal 
equilibrium and the Boltzmann equation (BE) to capture non-equilibrium phenomena. 

It is a major challenge to clarify the simplif actions and approximations posed in a micro- 
scopic theory which lead to such an effective model. We state spintronics as an example^'^ 
where spin-drift-diffusion equations proved to be a powerful tool for describing spin-polarized 
transport'^ and spin-transfer torques^'^ in magnetic mulitlayers. It has been demonstrated 
that these equations can be derived from a spinorial BE®"^. Hence, the missing link in a 
systematic, qualitative understanding is the derivation of a spinorial BE starting from a full 
quantum mechanical treatment. This is the main goal of the present study. 

Let us briefly discuss some well-established results in order to position the present work 
in an appropriate context: The generic form of the scalar BE is 

dtf-{hJh,r, = C{f), (1) 

where f{x,r],t), / > 0, is a probability distribution on the 2(i-dimensional phase space 
M^ X R^, h{x,T]) stands for the energy of a non-interacting particle and {/i, /}x-,r? = ^xh ■ 
Vj,/ — Vrjh ■ V xf denotes the Poisson bracket with respect to the position coordinate x and 
the momentum coordinate rj. The collision operator C on the right-hand-side (rhs) of (1) 
models short range interactions between particles or with obstacles, e.g. impurity centers or 
phonons in case of electronic transport in semiconductor devices. C is usually an integral 
operator and, moreover, non-linear in case that it describes interactions between identical 
particles or accounts for quantum statistics. Eq. (1) is referred to as the semiclassical BE 
since microscopic properties like the electronic bandstructure and quantum scattering rates 
can be described in terms of h{x,ri) and C{f), respectively. 

The incorporation of further quantum phenomena like coherence and entanglement^"^^ 
creates a need for either quantum corrections to the BE or quantum versions thereof, called 
quantum Boltzmann equations^^"^^. Moreover, the recent emergence of spintronics^^ raised 



the question of how to describe scattering of spin-coherent electron states in magnetic mul- 
tilayers or domain walls by means of a kinetic equation^' ^^"^^. In the spin-coherent regime, 
the BE (1) is replaced by 

dtF - {hi, F},,, + t[Q, F] = Q{F) , (2) 

where F{x,r],t), Q{x,ri,t) are hermitian 2x2 matrices defined on the phase space and 
[Q, F] = QF — FQ denotes the commutator. F is the distribution matrix, the eigenvalues 
of which give the scalar distribution functions of the two spin species. The term Q is an 
exchange field that mixes the two spin distributions. Equation (2) is referred to as the 
spinorial or matrix Boltzmann equation (SBE)^'^. Possanner and Negulescu^ studied linear 
collision operators Q which feature spin-dependent scattering rates, for example 



jdv' (^S''/'F{v')S''/' - ^SFiv) - If{v)S 



Q{F){^) = drj'i S"/'Firj')S''/' - -SF(r/) - -F(r/)5 , (3) 



Here, S = S{ri, rj') is a strictly positive, hermitian 2x2 matrix, whose eigenvalues denote 
the scattering rates from rj to t]' for the two spin species^^'^^ and 5" = Sit]',!]). The left- 
hand-side (Ihs) of Eq. (2) has been derived on a rigorous basis by Hajj^. A derivation of 
(3), which is able to relate the scattering matrices S* to a microscopic Hamiltonian will be 
accomplished in the course of this work. 

The derivation of QBEs or the SBE starts at the microscopic level by defining a suitable 
model Hamiltonian. Then, the natural framework to pass from the quantum to the kinetic 
level is the Wigner-Weyl formalism of quantum mechanics^^'^^. There exists a plethora of 
results regarding this passage for the scalar (spin-less) case, some of them we shall briefly 
mention here (for further information the reader is urged to view the references in the articles 
cited below). We remark that for the case that the eigenvalues of S are identical, performing 
the trace in Eq. (2) leads the scalar BE (1). QBEs have been obtained in the framework 
of generalized Kadanoff-Baym non-equilibrium Green's functions'^^'^'^'^^ and by a monitoring 
technique^^. On the rigorous level, the linear BE has been obtained from the single particle 
Schrodinger equation with a Gaussian random potential in the weak-coupling^^'^^ and in 
the low-density limit^^, respectively. The non-linear BE was derived by starting from the 
many-body Schrodinger equation with weak pair interaction potential and by studying the 
quantum version of the BBGKY-hierarchy.^^ 

In this work we apply a different strategy for passing to the kinetic level. Our starting 
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FIG. 1. (Color online) Schematic illustration of the strategies to pass from the quantum level to 
a semiclassical BE. The path employed in the present work is indicated by solid arrows while an 
alternative approach is indicated by dotted arrows. 

point for the semiclassical analysis will be a master equation of the Lindblad form^^ de- 
scribing a single quantum particle in contact with its environment^*^'^^. Semiclassical limits 
of Lindblad type master equations have already been considered. '^^ In particular, we shall 
start from the hierarchy of master equations derived by Possanner and Stickler^^. Master 
equations of Lindblad form describe the quantum evolution in terms of a semigroup law 
(quantum dynamical semigroups)^^, just as the BE does for the classical evolution on the 
kinetic level. Therefore, by starting the semiclassical analysis from the Lindblad equation 
instead of the von Neumann equation, the passage from the quantum to the kinetic level 
has been decomposed into two stages as sketched in Fig. 1: 

(1) In the quantum regime, one performs a Markovian limit that leads to dynamics de- 
scribed in terms of a quantum dynamical semigroup (Born-Markov limit^^), 

(2) in the Markovian regime, one performs the semiclassical limit (scaled h —^ 0) in order 
to obtain the BE. 

Corrections to the BE arise at each of the two stages. One obtains non-Markovian corrections 
at the first stage and quantum corrections in ascending powers of h (scaled) at the second 
stage. This paper deals solely with the second stage, while the first has been accomplished 
by Possanner and Stickler^^. 

This paper is organized as follows: in Sec. II we shall agree on some notations and 
specify the physical system under investigation. In Sec. Ill we explicitly evaluate the 
integral kernel of the dissipator, in Sec. IV we introduce the Wigner transform of the state 
operator and transform the whole master equation into the Wigner representation. Here we 



derive a Wigner equation equipped with a linear collision operator which features momentum 
relaxation as well as spin decoherence. The details of the derivation are explicated in App. 
C - App. E. In Sec. V we discuss the quantum collision operator and, finally, in Sec. VI we 
introduce the semiclassical scaling and define the different semiclassical scenarios which we 
will regard in this work. Moreover, we draw the semiclassical limit for these scenarios and, 
thus, derive the spinorial BE with a collision operator of the form (3) as well as the Bloch 
equations. Conclusions are drawn in Sec. VII. 

II. NOTATIONS AND MODELING 

We consider the evolution of a single quantum particle with two spin degrees of free- 
dom (spin 1/2-particle), henceforth called the 'system'. This particle interacts with the 
'environment' B which in turn is composed of A^ identical spin 1/2-particles. 

The dynamics of the system's particle are governed by a master equation of the Lindblad 
form^'^, 

d,p = -'-[Kf,p]+V{p). (4) 



Here, p is the density matrix defined on a one-particle Hilbert space "H, h denotes the Planck 

Vm 

— '-'— ^t^ - acxnig m tt reaas iiq ■ = uq -\- rimf, 



constant, -f^o and V stand for the system Hamiltonian and the dissipator, respectively. The 
Hamilton operator -ff^ acting in "H reads H^ = Hq + Hmj, where Hq is a one-particle 



Hamiltonian and one defines the mean-field operator 

Hmf = tiB (HiI <^xb) . (5) 



Here, Hi is an operator acting in the composite Hilbert space l-i ® I-Lbi which describes 
the interaction between the system's particle and the environment (or bath), 1 acting on 
T-i is the unity operator in l-L and xb is a predefined equilibrium density matrix on T-Lb- 
The operation tr^ (■) stands for taking the trace over the degrees of freedom of the bath. 
Moreover, the action of the dissipator T) is defined by 

Zif,,^ Ah-™/ lir'^f 



^(■):=-^tr^([^rJ^r,-®XB]]j, (6) 

where tq denotes the characteristic timescale'^^ of the system's dynamics -f^^ and H^ = 
Hj — Hmf ® Ifi- Please note that in writing (6) we assume that the Hamiltonian ffj is 
associated with the same characteristic energy eo as H^, i.e. we rescaled Hj"^ — )■ H^^eo/ej 



where ej denotes the characteristic mean- field corrected interaction energ^^. Equations (4)- 
(6) are vahd for very fast relaxation of bath states towards xb and ti/tq ^ 1. It has to 
be emphasized that due to the assumptions incorporated in the derivation of Eq. (4) we 
restrict our discussion to a case in which the system's particle is distinguishable from the 
particles constituting the environment. Furthermore, we note that under certain premises, 
Eq. (4) may account for the dynamics of the system's particle towards a unique equilibrium 
state.36 

Let us briefiy comment on the physical picture employed: We assume that Eq. (4) 
provides a proper description of the quantum dynamics of the system's particle in contact 
with its environment. It is the aim of this work to draw the semiclassical limit of Eq. (4), 
i.e. to regard the dynamics of the system's particle in a regime in which quantum effects 
cease to be observable. This goal is achieved in three steps: in a first step we shall rewrite 
Eq. (4) as an equation for the integral kernel p{x,x',t) = {x \p(t) \x') in position space, in 
a second step we shall derive the Wigner representation of Eq. (4) and, finally, in a third 
step we shall draw the semiclassical limit of Eq. (4) . 

However, we need to clarify some notations first. In what follows we shall denote the 
position and momentum coordinates of the system's particle by x G M^ and r] G M^, respec- 
tively, and the position coordinate of the n-th particle in the environment -B by 2„ G Mf . 
Moreover, we introduce the short-hand notation Z = {zi,Z2, ■ ■ ■ ,zn) G M.^'^. For multiple 
integrals we use the abbreviations 



dZ = / d^i / dz2 • • • / dzN , (7a) 

dZn = dzi... / dzn-i / dzn+i ... I dz^ , (7b) 



/ dZnm = dZi... dZn-1 / dZn+1 ... i dZm^i j dZm+1 • • • / d^Ar . (7c) 



The spin degrees of freedom of the system will be labeled by roman lower case letters, such 
as z,j G (1,2), while the spin degrees of freedom of the n-th particle in B are labeled by 
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Greek letters, i.e. an,Pn ^ (I72). In what follows we shall use the multi-index notations 

{a} = (ai,a2,- • • , "at) , (8a) 

{«}„ = («!, . . . , a„_i, an+i, . . . , Oat) , (8b) 

{0i}nm = («1, • • • , ttn-l, "n+l, • • • , "m-l, "m+l, • • • ^n) , (8c) 

{a, /3„} = (ai, . . . , a„_i, /3„, «„+!, . . . , a^v) . (8d) 

In particular, for the mult i- indices (8a)-(8d) we employ the abbreviation "^ta} — X]« « ■■■« 
and write the matrix elements of a 2^ x 2^ matrix T as T^°'^'^^^ = Ta-ia2-ai^,/3i/32---i3N- 

In what follows the set of hermitian 2x2 matrices is termed ^^(C) and we shall use the 
notation a = (o"i, (72, as), where 

^'^(Jo)' "^(°r)' "^(o-i)' *'' 

are the three Pauli matrices. Hence, any matrix G G J^{C) can be written in the Pauli 
basis (1,(?) with coefficients go eM. and g = {gi,g2,gz) £ IR^, respectively, 

G' = ^ol + ^-^, (10a) 

^o = ^tr(G), ^=^tr(aG'), (10b) 

where 1 is the 2x2 unity matrix. The eigenvalues of G can be expressed as (^o i 1^1 5 hence 
we call 1^1 the spin polarization of G. Moreover, we refer to g as the spin (spinorial) part 
or directional spin polarization and to g/\g\ as the direction of spin polarization of G.'^^ 
Moreover, we remark the following property: let A,B & ,^(C), then with the help of the 
decomposition (10) and the properties of the Pauli matrices we have 

i[A,B]=i[a-a,h-a] = 2{hxa)-a, (11) 

where a and h denote the spinorial part of A and B, respectively. 

With the help of these notations, let us further concretize the physics of the system 
under investigation. In the following we denote operators with a hat, such as A, and 
their integral kernels (matrix elements in the position space basis) without the hat, i.e. 
A{x^y) = (x A y). The integral kernel of Hq is assumed to be of the general form 
X Hn y) = Ho{x,y) G J^{C). Similarly, for the bath reference state xb we write 



viV 



{Z \xb\Z') = xb{Z,Z') G <S>n=i'^(^)^ ^■^- ^^ ^^ ^ hermitian 2^ x 2^ matrix at every 
point {Z, Z') G Mf^*^ x M^"^. The interaction Hamiltonian accounts for spin-dependent two 



particle interactions and is assumed to be diagonal in position space, ( x, Z 
Hi{x, Z, x', Z') = Hi{x, Z)6{x - x')6{Z - Z') G (g)^^ ^2(C), where 

N /n-l \ / N 

Hj{x, Z) = 5^ (g) 1 J ® V(a; - z„) ® (g) 1^ 



Hr 



X 5 Z; 



(12) 



n=l \fc=l / \fc=n+l 

Here, 1^ denotes the 2x2 identity matrix referring to the A;-th bath particle and V(r) G 
J^ ® J^2 stands for the spin-dependent pair interaction potential that depends on the 
distance r > between the system's particle and one bath particle as well as on their spins. 
With the help of the multi-index notation, we can write the matrix elements of Eq. (12) as 

AT 



Hj (x, Z) — 2_^ 5ai/3i ■ ■ ■ ^a„_i/3„_iV'Q,„/3„(a; " Zn)Sa„+if5„+i ' ' ' ^ajvfe) 

n=l 



:i3) 



where Sij denotes Kronecker's S and we note that V^„^„(r) are 2x2 matrices in the system's 



spin degrees of freedom which obey Vj 



an/3n ^ 



V8„a„(r) due to the hermiticity of V(r). 



The particular form of the pair interaction V(r) is given by the system of interest and, 
therefore, depends on the type of particles or quasiparticles constituting the system and the 
environment. Independent of its actual form we may express V(r) with the help of Eqs. 
(10) as 

V(r) = Voir) ® 1 + V{r) © a, (14) 

ai and we defined the hermitian matrices 
1 



where V{r) Q a = Y,i=i ^i( 



Voir) 
Viir) 
V2{r) 
Vsir) 



V^2ir) 



12 

VAr) 



(15a) 
(15b) 
(15c) 
(15d) 



[^ii(r) + V22(r)] 
Vi2{r) 

VUr) 
[Vu{r)-V22{r)]. 

Here, we already employed that V2i(r) = Vi2ir) and, hence, Vi{r) G J^{C). We shall 
frequently employ the representation (14) in what follows. 

Hence, we regard the dynamics of a single quantum particle in contact with its environ- 
ment. It is assumed that the dynamics are described in a proper fashion by the Lindblad Eq. 
(4) where the matrix element of the free particle Hammiltonian Hq is of the most general 
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form (x 



Hq y) = Ho{x^y) and the interaction Hj between the system's particle and the 
environment, Eq. (12), is composed of two particle interactions V(r) which is a function of 
the distance r between the interacting particles. 



III. INTEGRAL KERNEL OF THE DISSIPATOR 



It is the aim of this section to rewrite Eq. (4) as an equation for the integral kernel 
p{x,y,t) = {x \p{t) \y) G J^(C) in order to prepare the transformation of Eq. (4) into the 
Wigner representation in Sec. IV. Starting from Eq. (4), a straight- forward calculation 



gives 



dtp{x,y,t) + -C (h^^J {x,y,t) = Q{p){x,y,t), 



where we defined the integral kernel of the commutator 



C{A){x,y,t) = (x [A,p] 



y 



dx' [A{x, x')p{x', y, t) - p{x, x', t)A{x\ y)] 



(16) 



(17) 



X 



A 



x')G .^(C) and the integral kernel of the 



for some general operator A{x,x') = 
dissipator 

Qip){x,y,t) = {x\Vip)\y). (18) 

We note that Eq. (16) is entirely equivalent to Eq. (4) and will serve as basis to derive 
the Wigner equation in Sec. IV. The particular form of the dissipator Q as well as of the 
matrix element of the mean-field Hamiltonian Hmf{x) in Eq. (16) will be determined in 
what follows. For this sake, we insert into Eq. (6) the definition Hj = Hj — Hmj, see Sec. 
II, in order to obtain 



m ^ I 



tTB (2Hjp ® xbHi - HiHip ®xb-P® XbHiHj 

— '^HmfpHm.f + HmfHmfP + pHmfHm.f 



(19) 



To further simplify expression (19), we note that inserting the interaction Hamiltonian (12) 
gives contributions in which we can take the partial trace over all but one or two bath 
particles. It is therefore advantegeous to define the partial traces of xb over all but one or 
two bath particles, respectively, as 



{"}n 



(20a) 



X?i..„,J-n, z^) ■.= Y. j d^n^xlf '^™^^'""^(^, Z) , (20b) 

where the indices n, m still refer to specific particles. Since we assume that the bath particles 
are indistinguishable (Sec. II), we may omit these indices. It follows immediately that the 
corresponding matrices x i'^) ^ ^^(C) and x {z^^') ^ ^^(C) ® ^^(C) are normalized, 
i.e. 

Moreover, due to the indistinguishability it turns out to be beneficial to define the one- and 
two-particle spin-density matrices l>^^^\z) e .^(C) and N(2)(^,^') g J^(C) ® ^(C) with 
matrix elements 

n%{z):=Nx%{z), (21a) 

^S'/3/3'(^,^'):=iV'xi'i;3;3'(^,^'), (21b) 

respectively. Let us briefly comment on this definitions: Employing the notation (10), we 
call 

NW(2) =no(2)l+n(2)-<T, (22) 

the spinorial density of bath particles. Here the scalar part nQ{z) is the density of bath 
particles at z and the spin part n{z) is the directional spin polarization of the bath at z. 
The spin part n{z) is proportional to the magnetization m{z) of the environment'^^, thus, 
the spin-density matrix N^^)(z) at position z contains the complete spin resolved information 
about the probability of finding a bath particle at position 2. In a similar fashion we regard 
the matrix N^^^(2:, z') as the spin-resolved two-particle density matrix of the environment. 

With the help of the above definitions a straight forward calculation allows to express 
the mean-field interaction Hamiltonian Hm.f{x), 

Hr,^f{x) = f dzVo^^ix - z)n^^%) , (23) 

as well as the integral kernel of the dissipator (see App. A) 

9t r r 

Q{p){x,y,t) =—^ / dzdz'Kisfjfaa'{z,z') Vap{x - z)p{x,y,t)Va'p'{y ~ z') 

- 2^"/3(^ - z)Vo,^p'{x - z')p{x, y, t) - -p{x, y, t)Vafs{y - z)Va>p>{y - z') , 

(24) 
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for the system under investigation. Here and in what follows we shall employ Einstein's sum 
convention. Moreover, naa'i3(5'{z, z') are the matrix elements of the modified density-density 
covariance matrix K.{z, z') defined via 

K{z,z') = C{z,z') + -D{z,z'), (25) 

where D(2;, z') G J^{C) ® ^(C) is given via daa'(3i3'{z, z') := n^^{z)6aa'Si3(3i6{z — z') and 
we introduced the spinorial density-density covariance matrix (or environmental covariance 
matrix) 

C(z, z') = N(2)(^, z') - NW(2) ® N«(2')- (26) 

Hence, we determined all components of the evolution equation of the integral kernel 
p{x, y, t), Eq. (16). The particular form of the integral kernel of the commutator as well as 
of the dissipator were determined for the interaction potential (12) to be of the forms (17) and 
(24). The action of the dissipator Q on p is, thus, determined by the matrix elements of the 
modified spin-resolved density-density covariance matrix K(z, z') (25) and the interaction 
potential V(r). Furthermore, the explicit form of the mean-field Hamiltonian Hmf{x) was 
obtained, Eq. (23). The derivation of the Wigner equation in Sec. IV will essentially be 
based on these equations. 

However, before proceeding to the next section let us briefly discuss the integral kernel 
of the mean-field interaction, Hmf{x), Eq. (23), in more detail. The mean-field interaction 
Hmf{x) can be regarded as the partial trace tri(-) over the bath's particle spin degrees of 
freedom of the convolution of the potential V(r), Eq. (12), with the spinorial density N^^^ (z), 

Hmfix) = fdztii [Y{x - z)NW(^)] 

dzV,{x-z)ni{z), (27) 

where the sum runs from i = 0, . . . ,3 and the matrices Vi{r) G J^(C) have been defined in 
Eq. (15). 

The mean field interaction (27) is easily concretized for two particularly interesting physi- 
cal situations. As a first example we regard the case that the interaction is spin independent, 
i.e. 

Y{r) = v{r)l®l, (28) 

where v{r) ^M. and we obtain the familiar expression 

H^,f{x) = / dzv{x - z)no{z)l, (29) 
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which is known from spin independent mean- field theory^^. Here the mean- field interaction 
is substantially determined by the total density of bath particles no{z) at position z and, 
therefore, insensible for any spin-polarization of the environment. On the other hand, in the 
case that 

V(r) = t;(r)(T3 ® ^3, (30) 

i.e. a distance-dependent 3-polarized spin-spin interaction, we have 

Hmf{x) = / dzv{x - z)ns{z)a3, (31) 

which is a mean-field interaction that depends on the 3-spin polarization of the environment 
n^i^z) at position z. Since the spin polarization of the bath may be connected to a magne- 
tization m^lzY^ via m^lz) oc n^lz), we interpret Eq. (31) as the interaction of the system's 
spin with a position dependent effective magnetic field in 3-direction. Since we may write 
the mean-field interaction as Hmf{x) = h!^ {x)a-i we regard /i™ (x) as the 3-component of 
the effective magnetic field at position x. Here, the 3-component of the effective magnetic 
field /13 (x) is given by a convolution between v{x — z) and the magnetization m^z) oc n^{z). 
We shall come back to these two particular examples in the course of the following 
sections. 

IV. WIGNER REPRESENTATION OF THE LINDBLAD EQUATION 

It is the aim of this section to apply the Wigner transform W to Eq. (16). The ensuing 
transport equation for W = W(p) is a spinorial Wigner equation equipped with a quantum 
collision operator stemming from the dissipator Q defined in Eq. (24). This result will serve 
as a starting point for the semiclassical analysis performed in Sec. VI. 

In what follows we shall refrain from an explicit notation of the time argument t. Let us 
define the Wigner transform of p{x,y) G .^(C) and its inverse. The element- wise Wigner 
transform reads 

Wij{x,r]) := W{p,j) = — — ^ J dypij (^ + f'^^ " |) exp l-^y ■ vj , (32) 



and the corresponding inverse transform is given by 

/ / T* I 7/ 

Pij{x,y) = W~^{w^j) = / dr]Wij -^T—,V ) exp 



-(x-y) -r] 



(33) 
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By convention, we shall understand all Wigner transforms in an element-wise fashion as 
defined in Eq. (32) and (33), i.e. we denote W = yV{p). We remark that the Wigner function 
W must not be interpreted as a phase space distribution function since its eigenvalues also 
take on negative values. The size of negative regions of the Wigner function display the 
wavefunction's ability to interfere and may be used as a measure for the non-classicality of 
the system. ^° 

We now apply the Wigner transform (32) to Eq. (16) for the integral kernel p. For this 
sake, let a(x, rj) = {27ihYW[A{x, y)] G ^(C) denote the phase space symbol of an operator 
A acting in the system's Hilbert space "H. Then we write for any phase space symbol a the 
Moyal bracket^^ as (see App. B) 



Cn{a){x,7],t) ■.= W[C{A)]{x,r],t) 

lp|d.d.'|dede'[a(/ + ^.,e + ^r 



{2TThy 



a I z- --z,^--^' 



X exp 



4 



h 



z-{r]-C) 



exp 



h 



e' ■ (x - z') 



(34) 



which is the Wigner transform of the integral kernel of the commutator, C{A) defined in 
Eq. (17). In the general case, we remember that H^^ {x,y) = HQ{x,y) + Hmf{x,y), where 
the mean-field interaction Hmf{x,y) = Hmf{x)6{x — y) was concretized in Eq. (23). Let 
^o{x, rf) = {2TrhyyV{Ho) denote the symbol of Hq in phase space and let f)Q = f)o + Hmf{x) 
be the phase space symbol of H^l^. Then the Wigner-transformed Eq. (16) may be written 
in a compact form as 



dtW+'-Cnii)o^] 



To 
^2 



Qh{W). 



(35) 
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where the Wigner representation of the dissipator (24) reads 

Qn{W) := W [Q{p)] 

= ^ / dx'd^d^' / dri'KfJi3>aa'{z,z') 



{27rhy 



X 



1 



Va^{x+-X' -z] W'Va'^' ( X 



-X 



-K 



a/3 



X+ -X' - z] Va'p' ix+-x'-z']W' 



-\w'V^p (x-]^x'-z] V^,p. ( X 



-X 



X exp 



--X' . iv 



(36) 



with W = W{x,ri',t). 

Eq. (35) is the Wigner equation for a single spin-1/2 particle in interaction with its 
environment and is entirely equivalent to Eq. (16) of operator symbols in position space, 
from which it was derived. The Ihs of Eq. (35) represents the free flight of this particle 
and is therefore analogous to the equation obtained from the von Neumann equation in the 
Wigner picture, except for the mean-field correction. For instance, if H^ (x, y) gives the 
symbol [)o (x,?]) = [?7^/(2m) + u{x)] 1 + Q,{x,rj) ■ a + H^f{x), with the particle's mass m, 
the scalar external potential u{x) and the exchange field Q{x,ri) ■ a, one obtains the free 
flight term of the Vlasov equation 



n 



A ^0 



,"i/ 



m 



■ V^W + ^Cn{ul) + ^Cn [n ■ a + H^f) . 



(37) 



Here, Cfi [ul) is of the well known form (see App. B) 



U{ul) 



{2nhY 



X exp 



dx' / drj' 



u \ X + 



X 



u X 



X 

~2 



w 



--X' . (, 



V 



(38) 



The mean-field correction in Eq. (37) is a first contribution stemming from the interaction 
between the system's particle and the environment and its particular action depends on 
the form of the two particle interaction V(r). As illustrated in Eqs. (29) and (31) the 
mean-field interaction may give rise to a scalar as well as to a spinorial contribution to the 
one-particle Hamiltonian HQ{x,y). In particular, if V(r) is spin-independent, i.e. of the 
form Eq. (28), the mean field contribution to Eq. (37) is of the form Eq. (38) and H.mf{x) 
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behaves like an external scalar potential u{x). On the other hand, if V(r) models a spin- 
dependent interaction of the form Eq. (30) a brief calculation demonstrates that the mean 
field contribution to Eq. (37) takes on the form 



^h ( hT-^as 



'W^f (^'\ 



drj'h^ {rj' 



cT;,W(x,7]-'^,t\ -H^U,r/ + |,t)a3 



X exp —X ■ r] 
n 



(39) 



Here, h^ (t]') denotes the Fourier transform of the 3-component of the vector hmf{z), see 
App. C This term may be interpreted in the sense that it is accounting for the mixing of 



momentum components of different spin species due to the effective magnetic field /i™ (z) 
induced by the spin degrees of freedom of the environment (Sec. III). 

The rhs of Eq. (35) is a quantum mechanical collision integral (or collision operator) for 
the Wigner equation. We note that Eq. (35) together with Eqs. (34) and (36) will serve as a 
starting point for the semiclassical analysis carried out in Sec. VI, however, let us comment 
on some general properties of the quantum collision operator first. 



V. THE QUANTUM COLLISION OPERATOR 

We discuss some general features of the collision operator Qh{W) acting on the spinorial 
Wigner function W according to Eq. (36) in order to understand its action in more detail. 
This will prove to be crucial in the semiclassical analysis of Eqs. (35), (34) and (36) in Sec. 

VI. We insert inverse Fourier transforms (see App. C) of the interaction potential. 



Vafsix) = I d^V„^(0 exp {t^-^ 



(40a) 



and of the modified spin resolved density- density covariance. 



I^l3f}'aa'{z,z') = / dC,d^' KfJi3>aa'{C, O ^W 



^{z-^+z'-a 



(40b) 
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into Eq. (36) and subsequently integrate over x', z and z' in order to obtain 
Q,,(P^)(x,r/,t) = 2{2Txhf^ I ^i^^kf,fi,^^\i,0 



V^fi{i)W\x,^ 



t v^>p, ia 






X exp 



h 



x-{^ + 



(41) 



We denote this as the representation of the colhsion operator (36) as a momentum space 
integral Whether or not this representation is more convenient than Eq. (36) depends 
on the system under investigation. However, we shall now discuss the special case that 
K.{z, z') = K(2 — z') and it will turn out that in this case representation (41) is the more 
convenient one. We stress that it has to be checked carefully whether or not the simplification 
K(z, z') = K(z — z') is valid for the system under investigation. For the sake of a simplified 
discussion, we shall assume it to be a valid approximation in what follows. 

From definition (25) we remember that K.{z,z') = C{z,z') + Y){z,z'). For translational 
invariant systems the covariance C{z, z') = C{z — z'). If we further restrict our discussion 
to the case of a constant particle distribution, i.e. N'^^)(z) = const, we also have T){z,z') = 
Y){z — z'). It is, therefore, a sufficient condition to assume a space homogeneous environment 
whose spinorial density is constant in space in order to justify the above simplification. 
Hence, we regard the case that K(,^,,^') = K (2") ^(^ '^ O^ see App. C, and, therefore, 
Eq. (41) simplifies to 






(42) 



with W = W(x,r],t) and W = W(x,r]',t). Please note that we employed that Va/si^) = 
Vai3{—C,) according to Eq. (12). In what follows we shall highlight some general features of 
the quantum collision operator (42). 
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It is demonstrated in App. D that we may rewrite Eq. (42) in the even more convenient 
form 



Qn{W){x, r/, t) = 2{2nhf^ J dr]'p,{r] - r/' 
S,^ (r/ - r/O W'S, iv - 
-S, [r] - V) S, iv - v') W 



X [S, [r] - r/O W'S, (r/ - r]') 
1 
2 



(43) 

i.e. as a sum of coUision operators of the form of the Bohzmann equation, Eq. (3). Hence, 
in the particular case that the modified covariance fulfills K(2, z') = K(2 — z'), the quantum 
collision operator (36) turns out to be of the form proposed for the semiclassical SBE, Eq. 
(3). However, since W can take negative eigenvalues it is not a collision integral in the 
classical sense as already emphasized at the beginning of this section. Here, Pii^t]') G M, 
Si{ri') G J^(C) and z = 0, . . . ,3. As demonstrated in App. D the hermitian matrices Si 
are linear combinations of the matrices Vi, where the weights are determined by the matrix 
K(?7'). The scalar functions Pi{r]') G M are the eigenvalues of K (?]'). 

Let us briefly clarify the terminology used in the subsequent analysis: According to Eq. 
(10) we may express W a.sW = wl+tv-a, where we denote hjw = 1/2 tr (W) the scalar part 
and by uJ = 1/2 tr (Wa) the spin part of W. Furthermore, we call N = J drjW the spinorial 
particle density (or distribution) and M = J dxW the spinorial momentum distribution, 
which may also be decomposed according to Eq. (10). Again, the corresponding terms are 
referred to as the scalar and the spin part of N and M , respectively. 

The collision operator (43) may be decomposed according to 

Qn{-) = Q'i\-) + Q^'\-), (44) 

where 

Q'i\W) = 2{27rhy'' fdv'p^ir] - 7]')S, (r/ - r]') [W - W] S, (r/ - r/') . (45) 

and 

Q?\W){X, r/, t) = {2nhY'' J dVp.(^') ^ (V) ,W],S, (V)] • (46) 

We note that QJ^ is of the general master equation form ^gain term — loss term\ It is easily 



-h 
observed that we have 



JdvQi!\w)=(}, (47) 
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for arbitrary W, i.e. the spinorial particle density A^ is conserved by Qf^ ■ To be more 
specific, the particle distribution N{x, t) is not affected by the action of Q)^ while the scalar 
as well as the spin part of the momentum distribution M{ri, t) are changed. 

(2) 

In a similar fashion we obtain for Q)^ that 



tr 



)(2), 



QriW) =0, (48) 



(2) 



for arbitrary W, i.e. the scalar part w of the Wigner function W is conserved by Q\ . Thus 

(2) -» 

Q^ acts solely on the spin part w. Interestingly, the conservation of w signifies that both, 
the scalar particle distribution n{x, t) as well as the scalar momentum distribution ?7i(?7, t) 

(2) 

remain unaffacted. The operator Qj^ may therefore be identified as accounting solely for 
local spin-flip processes. 

Combining Eq. (48) and Eq. (47) gives the property 

fd7]tr[Qn{W)] = 0, (49) 

for arbitrary W. Thus the scalar part of the spinorial particle density, n{x,t) is conserved 
by Qh as in the well-known scalar case.^^ However, the spin part and in particular the local 
spin polarization \n{x,t)\ is changed, i.e. Qn accounts for spin-flip processes due to its part 

We now investigate in more detail the operator QJ^ . Due to its particular form (46) 
we may regard the arguments x, rj and t of VT as parameters and assume that the domain 

(2) 

of Q)j is J^{C) in the following analysis. At first, we are interested in its kernel denoted 
by Ker(Qj-j ), because, as will become apparant in what follows, we may under certain 
premises reason from the structure of Ker(Q^ ) on the stationary spin distribution of W for 
some given V(r). The Ker(Q|j ) is defined as the set of all matrices A G ^^(C) for which 
Q^ (A) = 0. As demonstrated in App. E we have 

Qf\A) = 0^[Sii7]'),A] = 0, \/t and Vr/', (50) 

if detailed balance is requred for A. Since Q}^ is linear and self-adjoint with respect to 

(2^ 

the Hilbert-Schmitdt scalar product ti (AB), A,B& .^(C), its domain D{Q\ ) may be 
decomposed into 

^(Qf) = =^2(C) = Ker(Qf ) © Ker(Qf )^ (51) 
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where 

Ker(Qf )^ := |a G ^2(C)|tr (AE) = 0V5 G Ker(Qf )} , (52) 

is the space orthogonal to Ker(Q^ ). 

If we knew the projection of some particular state W on Ker(Q)j ), we could assume that 

(2) 

Qj^ can be approximated by a relaxation time ansatz 

_(2) W -w 

Q'i\w):='-^-^^, (53) 

where W denotes the projection of W on the kernel of Q\ for some given V(r) and T2 
is the mean relaxation time. The action of the operator Q^^ as defined in Eq. (53) on 
some function W is that it relaxes W into W, where tr (W) = tr (W^ since tr ( Q^ j =0. 
We stress that, again, whether or not Eq. (53) is a valid approximation has to be checked 
carefully for the problem of interest. 

It is the aim of the following paragraphs to determine the equilibrium distribution W G 
Ker ( QJj ) for a given two particle interaction potential V(r). First of all, we deduce from 
Eq. (50) and (11) that for A G Ker(Q^^^) we have 

[S^{r]'),A] = 2t[si{7]')xd]-a = 0. (54) 

Since A does not depend on 77', this is for Si{r]') ^ and a ^ only possible if the direction 
of the spin part of Si^f]') does not depenend on r]', i.e. Si{ri')/\si{ri')\ = const. On the other 
hand, if it does depend on r]', i.e. if Si{ri')/\si{ri')\ 7^ const, the only possible solution to Eq. 
(54) is a = 0. In particular, it follows that in this case the spin part of W is zero. Hence, 
the unique projection W \-^ W reads 

W=hT{W)l, (55) 

i.e. the operator Q^ destroys spin polarization on a finite time scale T2 whenever the 
interaction potential's direction of the spin polarization Si{'r]')/\si{ri')\ ^ const. In a similar 
fashion, we note that Eq. (55) is the sole solution to Eq. (50) if the directions of the 
spin parts of Si are not identical for different i. We highlight that this mechanism of spin 
depolarization is not based on spin orbit interaction^'^ but on direct spin-spin interaction 

between different particles. Again, we emphasize: Whenever the direction of the spin part 

(2) 

of the interaction Si is a function of r]' or not equivalent for different i, the action of Q^^ 

unavoidably leads to spin depolarization on a finite time scale, as described in Eq. (55). 
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We now concentrate on the special case that the direction of the spin part of Sii^r]') is 

(2) 

constant in r]' and identical for all i. If Si{ri') = then Q)-^ (■) = and the operator Qn does 
not account for spin-decoherence at all. However, if 3,(77') 7^ then Si{ri')/\si{ri')\ = const is 
identical for all i, i.e. we may write 

s.iv')=l^iv')\ (56) 

where 71(77') is some scalar function and |A| = 1. Hence, 

S,iv') = Siiv'n + l^iv')^ ■ <?, |A| = 1. (57) 

Let S be the matrix which diagonalizes X-a and, according to Eq. (50), also W G Ker(Q^ ). 
We then obtain that the unique equilibrium spin configuration W reads 

|y = Itr (VT) 1 + itr (asS^W^S) SasS^ (58) 

In order to prove the result Eq. (58) we note that we may decompose an arbitrary matrix 
W according to W^ = &{W) + ff{W), where ^{W) denotes the matrix whose off-diagonal 
elements are zero and ff{W) denotes the matrix whose diagonal elements are zero. Now, 
according to Eq. (50) we have 

W e Ker(Qf ^) ^ S^WS = ^(S^WS), (59) 

is solely diagonal. In a similar fashion 

W e Ker(QfV ^ S^W^'S = ^(S^VT'S), (60) 

is solely off- diagonal. ^^ Hence, with the help of Eq. (51), since S^IVS = ^(S^iyS) + 

^(S^VrS), we note that 

S^WS = ^(S^l^E), (61) 

defines the unique projection onto Ker(Q^ '). Finally, inserting into (61) the expression 

^(Et^yS) = -tr (W) 1 + -tr (CTaS^l^S) ^3, (62) 

and solving for W proves statement (58). 

Hence, it was possible to identify the equilibrium spin configuration W for a given two 

particle interaction potential V(r), provided that Qn, leads to a unique equilibrium distri- 

(2) 
bution in the kernel of Q^ ■ This analysis of the quantum collision operator will show to 
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be a crucial ingredient in the following section in order to derive the Bloch as well as the 
Boltzmann equation. Furthermore, in the above discussion we identified the necessary cri- 
teria which allow for the treatment of spin depolarization via a relaxation time operator of 
the form Eq. (53) together with W given by Eq. (55). Such an operator was employed by 
Possanner et al.^ to derive spin drift-diffusion equations from a semiclassical BE. 

However, before we proceed to the semiclassical analysis of the Wigner equation (35) let 
us briefly discuss the two particular examples already illustrated in Sees. Ill and IV. 

First, we regard an interaction potential of the form (28), i.e. a spin-independent inter- 
action. In this case Q\ vanishes for all W , i.e. spin polarization is not destroyed. Hence, 
the quantum collision operator (45) takes on the form 

Qt,{W) = jdv'uiir] - 7]')\v{7] - r]')\' [W{x,v',t) - W{x,v,t)] , (63) 

where we deflned 

a;i(r/') = 2$^p.(r/'). (64) 

i 

We note that the collision integral Eq. (63) independently acts on both spin directions and 
for each spin species it is of the classical scalar form, however, it has to be kept in mind that 
W is not a proper distribution function. The weight uJi{ri') is given by the matrix trace of 
the corrected density- density covariance K(?7') (note that the trace is basis independent). 

In a similar fashion we obtain for a spin-spin interaction with a potential of the form (30) 
the collision integral 

Qn{W) = j dr/'a;2(r/ - i) Viv - i)f [a^W{x, r/', t)a^ - W{x, r/, t)] , (65) 

where 

a;2(V) = 2(-l)°+"'/^a„w(V)- (66) 

In particular, according to Eq. (58) the projection of W on the kernel of Q\ is of the form 

W = ^tr {W) 1 + itr {a-iW) as, (67) 

i.e. all particles are spin-polarized in 3-direction in the equilibrium state. 

VI. SEMICLASSICAL LIMIT: SPINORIAL BOLTZMANN EQUATION 

We shall now study the Lindblad Equation (4) [or, equivalently, Eq. (16) or the Wigner 
equation (35)] in a regime in which quantum effects other then spin coherence cease to 
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be observable. The passage from the quantum (or microscopic) world to the classical (or 
macroscopic) world will be described in terms of a small parameter e, referred to as the 
semi classical parameter or scaled Planck constant. This parameter tends towards zero as 
quantum effects become more and more negligible. 

Before discussing the mathematical subtleties of this transition we shall connect the 
semiclassical limit to the following physical picture: suppose Aq is a characteristic length on 
which quantum effects dominate the physics of the system and let Vq = Aq be the volume of 
a thus defined quantum box Bq- Moreover, suppose Eq. (35) provides a proper description of 
the system in the quantum box Bq. In the semiclassical picture we are interested in solving 
Eq. (35) in a domain Bs with volume Vs = Vq/e'^ for 5 <^ 1 with suitable intial and boundary 
conditions. Hence, e~'^ is a measure of how many quantum boxes with volume Vq fit into 
the regarded domain Bs- All wavefunctions are normalized within the domain Bs- The limit 
e — 7- can be understood as zooming out of the microscopic world Sq, such that wave- 
characteristics of particles are on a very small length scale compared to the macroscopic 
domain Bs, i.e. Vq/Vs = e'^ — )■ 0. 

On the other hand, we define a microscopic time scale tq on which quantum effects are 
dominant and which defines the corresponding quantum time domain To = [0,to]. Then, 
in the semiclassical limit, we regard a timescale r^ = tq/e characterizing the macroscopic 
time domain Ts- With the help of Aq and tq we may define the characteristic scales of the 
microscopic world i5o x To by the relations 

^ = ^ = 1, (68) 

n n 

where eo and ttq stand for the energy scale and the momentum scale associated with the 
microscopic world Bq x To, i-e. with Eq. (35). It is the aim of the following scaling 
considerations to express Eq. (35) in variables of the macroscopic world Bg x %. 

This is performed in three stages: (i) all functions f{x,ri,t) which appear in Eq. (35) 
are rewritten in such a form that we explicitly emphasize their characteristic wavelength 
and amplitudes. In particular, we rewrite f{x, rj, t) = Ocf if-,^,^), where ttc denotes the 
characteristic amplitude and Ac, vTc and Tc denote characteristic length-, momentum- and 
time scales. Thus, the function f{x,ri,t) is of order one with gradients of order one. 

In a second step, we introduce dimensionless variables 

t' = -, x' = f, r^' = ^ E' = ^, (69) 

Ts As TTg eg 
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where we introduced the semiclassical scales of the macroscopic world Bg x %■ According 
to the above considerations the semiclassical scales are connected to the quantum scales via 

^ = ^ = -• (70) 

As Tg 

It follows from the definition (68) of quantum scales that the above definition implies that 



(71) 



with Ss = eo and vr^ = ttq, i.e. the energy and momentum scales remain unaffacted while, by 
expressing Eq. (35) in the dimensionless variables (69), we regard a long-time, large-scale 
limit. 

In the third and final step, we define the characteristic lengthscale as well as the char- 
acteristic amplitude of all quantities which appear in Eq. (35) by posing suitable scaling 
assumptions. These assumptions determine the properties of the system under investigation. 
To be specific we assume in the following semiclassical analysis: 

i)l^\x,v) = eoi) (^A]t + eeon f f , -) ■ c?, (72a) 

yAs TTsJ yAs T^sJ 

V(^)(r) = aeoV(^), (72b) 

N<"(.) ^ ^N« (£) . (720 

where we omitted the overlines on the rescaled functions, however, indexed all quantities 
by (s) which implicitely depend on e. Further, a, b and i are scaling parameters which will 
assume different definite values in the following subsections. Depending on the particular 
choice of these scaling parameters we shall derive different macroscopic transport models. 

Let us briefly interpret the scaling assumptions (72). The free particle Hamiltonian f)o 
varies on the macroscopic scale A^. We note that the weak scaling of the spin part Q ■ a 
of the free particle's Hamiltonian {)o is necessary in order to preserve spin coherence in the 
limit e — !■ 0, as it was demonstrated by Hajj^. The particle distribution N*^^)(z) varies on the 
macroscopic scale A^. The scale on which the interaction potential V(r) varies is denoted by 
i. For instance, for I = Aq, Eqs. (72) describe short-range interactions between the system's 
particle and an environment whose density N'^^^(z) varies on the macroscopic scale A^. On 
the other hand, for i = Xs Eq. (72b) accounts for long-range interactions. The parameters a 
and b are the characteristic amplitudes of V(r) and N'^^)(2;), respectively. For instance, the 
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combination {a = 1, b = e) is referred to as the 'low-density scaling'^^'^^ while the scenario 
[a = e, b = 1) is denoted the 'weak-coupling scaling'^^'^^. 

In what follows we shall denote the rescaled Wigner distribution matrix W by W^^"^ 
in order to emphasize its dependence on the semiclassical parameter e. A short calculation 
based on the steps illustrated above shows that the rescaled Wigner equation (35) is obtained 
by replacing To/ff in front of the collision integral as well as 1/h in front of the commutator 
by 1/e. This follows from changing to the variables Eq. (69) and employing definition (71). 
In a similar fashion, all h-s appearing in Eqs. (34) and (36) are replaced by e. We obtain 
the rescaled version of Eq. (35) as 

where L^ and Q^ denote the rescaled operators (34) and (36) and \)^ = (Jq + -f^^/- ^'^^ 
e = 1, the domain on which Eq. (73) is defined is equivalent to the microscopic world i3o x To 
while for e — )■ we approach the macroscopic world Bg x Xs- 

We assume that a solution VT^^^ of Eq. (73), for suitable intial and boundary conditions, 
can be written as W^'^^ = F + eWi + . . ., where F{x,r],t) E ^(C) is a positive definite, 
normalized, slowly varying function of x, t] and t, called the spinorial distribution function. 
This assumption is necessary in order to obtain a well posed semiclassical transport equation 
because then for e — )■ 0, W^'^^ — )■ F is a proper distribution matrix in the classical sense. 

Before we refer to the specific limits, we shall briefly regard the system's particle free 
Hamiltonian [)q since this part is independent of the scaling parameters {a,b,i). It is 
demonstrated in App. F that one obtains 

Z y I- J x,r) I. J x,ri, 

+0(e), (74) 

which is commonly referred to as Moyal's bracket^^. From the scaling assumption (72a) 
we deduce that in zeroth order the spin part ^2 of [}q appears in the commutator while 
the scalar part f)l enters Poisson's bracket. Hence, provided the semiclassical limit of the 
mean-field term Cs and of the quantum collision operator Q^ exist, we may write Eq. (73) 
in the macroscopic world Bg ^ Ts, i.e. for e — > as 

d^F - {i)l, F}^^^ + tp ■a,F] + L{F) = Q{F), (75) 
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where L{F) denotes the seiniclassical hmit of the mean-field term and Q{F) the semiclassical 
hmit of the colhsion integral, see Eq. (73). Please note that the quantity F in Eq. (75) is 
a proper distribution matrix, i.e. Eq. (75) represents a macroscopic transport equation. It 
is the aim of the following subsections to specify the particular form of L(-) and Q(-) for 
different scaling scenarios. 

The requirement that L{F) must not diverge, i.e. that Eq. (75) exists, poses constraints 
on the different combinations of free scaling parameters (a, b, £). In particular, inserting 
Eqs. (72) together with the definition Eq. (69) into the rescaled mean-field interaction (23) 
gives 



^S(^) 



ab 



dznfl (z) V^,3 






(76) 



Let us briefiy discuss the two different scenarios for i for arbitrary a, b: (i) If ^ = A^ the 
above integral is of order ab. (ii) If £ = Aq the interaction F[^\{x) is strongly varying in 



X 



and N'^-'^)(2) enters only as constant, i.e. N(^)(0). We remark that the contribution of a 
strongly varying mean-field interaction H^Ax) to Eq. (75) vanishes in the semiclassical 
limit e — )■ for most cases. We shall come back to this point in Subsec. VIA. Due to 
the prefactor of 1/e in front of C^ \^m^ ) ^^ -^^- C''^) ^^ require that ab = e. Particularly 
interesting are the two specific cases (a = 1, 6 = e) and (a = e, 6 = 1), i.e. the low-density 
and the weak-coupling scaling. 



A. Short-range interactions 



Ao) 



Here, we study the semiclassical limit of Eq. (35) under the scaling assumption Eq. (72) 
for (a = 1, 6 = e, £ = Aq), i.e. a short-range low density scaling and (a = e, 6 = 1, 
i = Ao), i.e. a short-range weak coupling scaling. It is important to realize that we restrict 
our discussion to potentials which are integrable in M*^, i.e. Vapir) G L^(]R'^), such as the 
Yukawa potential (screened Coulomb potential). Please note that this excludes the bare 
Coulomb interaction. However, since in this case V{r/e) = eV{r), the strongly varying 
Coulomb interaction is equivalent to the slowly varying Coulomb interaction in the weak 
coupling scaling and will therefore be discussed in Subsec. VI B. 

First, we note from Eqs. (76) that the mean-field contribution vanishes in both scalinges, 
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i.e. L{F) = 0, since 

H^^lix) =^J<^^J d^4liez)Vo^piO exp(-2^ ■ exp (^x ■ e) ^ 0, (77) 

due to the Riemann-Lebesgue lemma and V^^(r) G L^(M.'^). Here, the Fourier transform 
Vai3{^) of Vai3{x/e) is independent of e, see App. C. 

In order to evaluate the semi classical limit of the collision integral, we make the following 
hypothesis on the environmental covariance C defined in Eq. (26), 



= 7r. {z - z') + ^ro ^^ j , (78) 

i.e. we regard a space-homogeneous environment whose covariance has a microscopic and a 
macroscopic part, Fq and F^, which are assumed to be behave as 7 when approaching the 
macroscopic regime. 

In the low density scaling, the resulting collision integral is for 7 = e of the form 

Q(F)=go(F), (79) 

where we have (App. G) 

Qo{F) = J dr/'pf ^(^ - r]') [S, (r/ - r/') F'S, {r^ - V) 
1 



^Si ij] - ri') Si {ji - ri') F 
lFS,{v-v')Sdv-v') 



(80) 



with F = F{x,r],t) and F' = F{x,ri',t). Here, the matrices Si{C,) are, again, hermitian 
linear combinations of the matrices Vi which are obtained by diagonalizing the matrix Fq, 
see App. D. The scalars p^ G M are the eigenvalues of Fq. In the low density scaling, Qq 
vanishes for 7 = e^ or weaker while it diverges for 7 = 1 or stronger. The same collision 
integral (80) is obtained in the short-range weak coupling scaling for 7 = 1/e. 

We note that the operator (80) is a linear collision integral of Boltzmann-form, i.e. in 
the case of a strongly varying density-density covariance C^"^-* we obtain the semiclassical 
Boltzmann equation. The complete discussion of the collision operator can be adapted from 
Sec. V, however, in the present case it was not necessary to assume that the bath density 
N*^^) = const. A spinorial Boltzmann equation of such a form has already been postulated 
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by Possanner et al.^, however, it has not been employed as a basis for the derivation of spin 
drift-diffusion equations. The coUision operator used to derive drift- diffusion equations can 
be obtained from (80) by rewriting the colhsion operator as sum of an operator accounting 
for momentum relaxation Qi and a second operator accounting for spin decoherence Q2, see 
Sec. V. The latter is then replaced by a relaxation time approximation of the form Eq. 
eqrefeq:relaxapprox where the projection of F onto the kernel of Q2 has to be assumed to 
be of the form F = |tr (F) 1. Within this work it was possible to identify the necessary 
criteria allowing for such a treatment, see Sec. V. 



B. Long-range interactions {£ = As) 

Let us regard an interaction potential varying on the macroscopic scale, i.e. i = Xg in 
Eq. (72). Again, we shall regard the low density {a = 1, b = e, i = A^) as well as the 
weak coupling {a = e, b = 1, i = As) case. Please note that the strongly varying Coulomb 
potential is equivalent to the slowly varying Coulomb potential in the weak coupling scaling. 
From Eq. (76) we deduce for both scalings 

L{F)=z[Hmf,F], (81) 

where H^f is the macroscopic mean- field interaction of the form (23), i.e. independent of e. 

In the case of long range interactions we employ the following ansatz for the environmental 
covariance C: 



C^'\z, z') = -fC^'\z, z') 



7 ( z z' 



e^^ \e e 



7rs(^,^') + 47ro -'- ' (82) 



i.e. in contrast to Eq. (78) we do not restrict to the space homogeneous case. 

In the low density scaling, the collision operator Q{F) takes on the form (App. H) 

Q{F) = Qld{F) = Qs{F) + g„(F), (83) 
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where for 7 = 5 



X VafS {X - Z) FVa'fi' {x - z') 
-2^Q/3 i^ - Z) Va'l3' {X - z!) F 
--FVap {X - Z) V^ipi {x - z!) 

stems from the slowly varying part of the environmental covariance and 

2 



(84a) 



Qn{F) 



(27r)^ 



dzn-ll{z) 



V^R (x - z) FVo^R ix - z) 



-2^a/3 i.^ - Z) ^«'/3' {X - Z)F 
-^FVap {X - Z) Va'p' {X - Z) 



(84b) 



stems from the matrix 0^"^^ In the long-range weak coupling limit Q{F) reads for 7 = 1/e 
(App. H) 

Q{F) = Qwc{F) = Qs{F). (85) 

Please note that we refrained from writing Eqs. (84a) and (84b) as momentum space in- 
tegrals and in the basis described in App. D for the sake of a more transparent notation. 
Of course, we may obtain such a representation by replacing the interaction potential by 
its (macroscopic) Fourier representation and carrying out the steps in App. D. (In par- 
ticular, for Eq. (84a) we note that in the case that the environmental covariance obeys 



Ts{z, z') = Ts{z — z') the representation as a momentum space integral will turn out to be 



more convenient. Moreover, the form of Qs{F) will in this case be very similar to Eq. (80), 
i.e. one obtains a Boltzmann collision integral.) 

Let us briefly discuss the ensuing transport equation in the long-range weak coupling 
limit for 7 = l/e. According to the above discussion we have 



dtF - {[)1, F}^ ,, + t[n-a + Hmf, F] = g,(F), 

The transport equation of the spin part / of F is of the form 

dtf + (V.f) ■ V,/ - V,{) ■ v./) + 2/ X (n + hmf 
= - J dzdz''y'i''\z, z')si{x, z, z') x Si{x, z, z') x / 



(86) 



(87) 
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where we employed Eqs. (10) and (11) and rewrote Qs with the help of the steps outlined 
in App. D. In particular, the vectors Si{x,z,z') are linear combinations of the spin parts 
Vi{x — z) of Vi{x — z), see Eq. (15), however, the weights may be functions of z and z' . The 
quantities 7j (;z, z') are the eigenvalues of Ts[z, z'). Equation (87) describes the precession 
of the directional spin polarization / under the influence of an external field Vt and an 
interaction with the environment. This interaction induces an additional field hmf and gives 
rise to the dissipator Qs-, which relaxes the vector / into a predefined direction, see Sec. V. 

Let us consider a spin located at a certain lattice point ( 'system ') which interacts with 
other spins located at different lattices sites ( 'environment^. Since the system's spin cannot 
move on the lattice the scalar part of the free particle Hamiltonian vanishes, [) = 0, i.e. 
it has no kinetic part and no scalar field is externally applied. We recognize that in this 
particular case Eq. (87) gives rise to the macroscopic Bloch equations of magnetism. To 
be more specific, we replace in Eq. (87) the dissipator Qs by a relaxation time ansatz as 
discussed in Sec. V and interpret the spin part / = l/2tr {Fa) as the magnetization of the 
system. 

We emphasize that it was, therefore, possible to derive the Boltzmann as well as the Bloch 
equations from the spinorial Wigner equation derived within this work. In a similar fashion 
one may impose further scaling assumptions and study the resulting classical transport 
models as well as their quantum corrections. 

VII. SUMMARY 

We derived a linear quantum collision operator for the spinorial Wigner equation. Fur- 
thermore, it was demonstrated that the Wigner equation gives rise to several linear semiclas- 
sical spin-transport models. We derived the Bloch equations as well as the linear Boltzmann 
equation as an example. Let us briefly summarize the main aspects of the derivation. 

We investigated the dissipative dynamics of a spin-1/2 quantum particle, referred to as the 
system, in contact with its environment, which is ni thermal equilibrium. It has been shown 
by Possanner and Stickler^^ that in the limit of vanishing system-environment correlations 
these dynamics are properly described by the Lindblad equation (4). The latter served as 
a basis of the current study. The Wigner representation of the Lindblad equation (4) is a 
spinorial Wigner equation (35) equipped with a quantum collision operator (41). It is then 
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demonstrated that the latter can be cast into the form of a Boltzmann colhsion integral, 
Eq. (42), provided that the spinorial density of the environment is constant and that the 
spinorial density-density covariance of the environment is space homogeneous. The 2x2 
hermitian scattering matrices Si are uniquely determined by the spin-dependent two-particle 
interaction potential V(r) and by the modified spinorial density-density covariance H{z, z') 
of the environment, defined in Eq. (25). The eigenvalues of the scattering matrices Si are 
the scattering rates for the two spin species. 

Moreover, the quantum collision operator (42) is composed of two qualitatively rather 
different parts. The first part changes solely the spinorial momentum density of the system 
while the second part accounts for local spin flip processes. Hence, the second part modifies 
the local spin polarization of the system in a fashion uniquely determined by the scattering 
matrices Si and the eigenvalues of the modified spinorial density-density covariance K(2;, z'). 
Furthermore, it is possible to identify clear criteria under which the interaction between the 
system's particle and the environment leads to spin decoherence or even spin depolarization 
in the long time limit. 

Finally, we performed a semiclassical analysis of the spinorial Wigner equation (35), i.e. 
we regarded the dynamics of the system's particle in a regime in which quantum effects 
other than spin-coherence cease to be observable. We restricted our discussion to the well 
established low-density and weak-coupling limits. In principle, several semiclassical evolu- 
tion equations for a positive definite, hermitian distribution matrix F can be obtained. As 
two particularly interesting examples we note the derivation of the Bloch equations for long 
range interactions and the spinorial Boltzmann equation for short range interactions and a 
spatially homogeneous environmental density-density covariance. ^^ This form of the Boltz- 
mann equation has already been used for deriving spin- coherent drift-diffusion equations in 
magnetic multilayers^. 

In summary, we remark that within this work it was possible to systematically establish 
the link between a full quantum-mechanical treatment of a composite quantum system by 
means of the von-Neumann equation and macroscopic linear spin-transport models featur- 
ing dissipation such as the spin drift-diffusion models. This makes the derived equations 
particularly interesting for applications involving graphene^^ (pseudo-spin formalism) or 
magnetically doped semiconductors^'^. 

possible to augment the resulting transport models with corrections which result from 
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the Born-Markov limit as well as from the semiclassical limit if deemed necessary. Within 
this work we restricted to the zeroth order equations in both scalings but the evaluation of 
higher order corrections is straight-forward. 
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Appendix A: Derivation of Eq. (24) 

We evaluate the matrix elements of the dissipator (19) by investigating the terms con- 
taining tr^ ( HjHjt (8) X-B ) • III what follows we employ Einstein's summation convention in 
order ot simplify the notation. We obtain 

tr^ (^H,HA ® xb) {x,y) = J dzdz' [n{N - l)xg]w(^, z') + Nxf^{z)6^^>6fsp'S{z - z') 

xVo,p{x-z)V^>p>{y- z')5{x-y) (Al) 

where Xa^'aa'^^^ ^') ^ ^ stems from Eq. (20b) for identical particles (indices n,m omitted). 
In a similar fashion one obtains the following relations 

trs (Hip ® xbHj^ {x,y) = N{N - I) j dzdz'xfl^^,{z, z!)V^f,{x - z)p{x, y)V^>i3'{y - z') 

+N f dzxf^{z)V^p{x - z)p{x,y)V^p{y - z), (A2a) 



tiB [H]p ® Xb) {x, y) = N{N - 1) / dzdz'xfl^^,{z, z')V^f,{x - z)V^,p,{x - z')p{x, y) 

+N I dzx^^l{z)V^p{x - z)V^p{x - z)pix,y), (A2b) 



tiB [p^xbHPj {x,y) = N{N - 1) j dzdz'xfp,o.Az^,z')p{x,y)Vap{y - z)Va'p'{y - z') 

+N j dzx^f^l{z)p{x,y)V^(s{y - z)V^p{y - z), (A2c) 

H„,fpHmf) (x, y) = N^ j dzdz'x^^l{z)x%,{z')V^p{x - z)p{x, y)V^'p,{y - z'), (A2d) 
Hlfp) (x, y) = N' j dzdz'x^^l{z)x^^Uz')Va^ix - z)V^,p:{x - z')p{x, y), (A2e) 
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and 



pHlA {x,y) = N' / dzdz'x^^i{z)xZ'i^')p{x,y)V^p{y~z)V^,piy-z'). 



(A2f) 



Using that A^(A^ — 1) f« A^^, applying definitions (21) together with (25) and (26) and 
inserting the relations (A2) into Eq. (19) yields the final result Eq. (24). 

Appendix B: Derivation of the Moyal bracket 



Within this appendix we derive the phase space symbol of the commutator Eq. (34), 

Ap 



i.e. the Moyal bracket. We calculate the Wigner transform of the element ( x 
Ap) (x, y) for some operator A as 



W 



Ap) {x,y) 



— — ^ / dx'dzA lx + -x', z]plz,x--x'] exp{-ix' ■ r]) 



xW I ^y^--x',^'jexp 



x + z I , . 






X exp 






exp ( ——X ■ 7] 



(Bl) 



where we used Eq. (33) and (x A y) = A{x,y) = (27r^) '^W ^[a{x,'r])]. We rearrange the 
exponential terms as 



^x-(e-n + ^^-(e'-o-^a:' 



^-2(^+0 



which suggests the substitution $,=$,—$,' and ^ = ^{C, + $,'). Hence, we have 

1 



W 



Ap) {x,y) 



{2TTh) 



2^ , dx'd2 / d^d^a ( ^^ + \x', e + ie 



xW { — -x,^--{)exp 



-(x-2;) ■^+^X' -{^-7]) 



.(B2) 



This expression may be rewritten in a more convenient form by replacing z' = ^^ and also 
^ = 2^. Finally, we obtain the first term of Eq. (34) as 



W 



Ap) {x,y) 



{2Txh) 



2d 



dx'dz' I didia{z' + ^x',i + ^i 



xW{z'-\x\i-\i]exV 



-{x- z') ■^+^x -{^-r]) 



■ (B3) 
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A similar calculation for the term ( x 

1 



pA 



y 



pA) {x,y) gives 



W 



pA (x, y) 



{2nh) 



2d 



dx'dz' I d^diw (z' + jx',^ + ji 



xa( z' - -x',^--^] exp 



-(x-z') ■^+^x' -{^-r]) 



(B4) 



Combining this relation with Eq. (B3) gives Eq. (34) as 



W[C(A)(x,y)](x,ri) 






a{z'-\z,i-\e 



exp 



-^' ■ (x - z') 



X exp 



-^-(^-0 



(B5) 



We shall now investigate the particular case (x 



B 

it follows from Eq. (32) that b{x,ri) = B{x). We retur to Eq. (Bl) in order to obtain 



y) = B{x,y) = B{x)6{x — y) where 



W 



Bp) {x,y) 



— — ^ / dx'B (x+-x'j p (x + -x',x--x' ) exp(-ix'-r/) 



and, therefore, for B[x) = 6(x)l the well known result Eq. (38) 

1 



j{i-v) -x' 



(B6) 



W[C{hl){x,y)] 



{2.hY ' ^^' ' ^^' 



6 ( X + -x' I — 6 I X — -x' 



ly' 



X exp 
follows, where we defined W = W{x,^',t). 



-^x'-(r/-0 



(B7) 



Appendix C: The Fourier transform 



We define the Fourier transform of an operator, which is diagonal in position space, i.e. 

V{x, y) = V{x)6{x — y) as 



and, therefore, its inverse as 



^|dxnx)exp(-lx.^ 



V{x) = I dr]V{r]) exp ( — x ■ rj 



(CI) 



(C2) 
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In a similar fashion we define the Fourier transform of a function K{x, y) which stems 
from a two particle operator i.e. K{x, y, x', y') = K{x, y)6{x — x')6{y — y') as 



m^o 



1 



dxdyK{x, y) exp 



In particular, if K{x, y) = K{x — y) we obtain 

1 



h 



[x-i + y-O 



(C3) 



I 



^^^' ^'^ = (2^ / '^'^^^''^ exp ( --X ■ e ) Ki + o 

Moreover, we will frequently employ the identities 

1 



and 



5{x) 



6{v) 



(2vr) 



{2ny 



driexp{iri ■ x) 



dxexp(— i?7 ■ x), 



(C4) 



(C5) 



(C6) 



where S{-) denotes Dirac's delta distribution. 

The semiclassical Fourier transform is obtained by replacing all appearing h-s in Eq. (CI) 
by e. Hence, we have 



t/(^)(ry) = — / dxV^^\x) exp ( — x ■ rj 

{2TTey J \ e 



(C7) 



and, therefore, also 

V^'\x) = f dr]V^'\r]) exp (-x-rij, (C8) 

for a single particle operator which is diagonal in position space. Here, the index (5) signifies 
that V may still be a function of e. In the particular case that V^^\x) = V{x/e) we obtain 
the important result that 

1 



V^'\ri) 



(27re)'^ 

is independent of e. In a similar fashion we obtain that 

1 



dxV ( — j exp I — X ■ T] ] = V{rj), 



i^(^)(e,0 



(2vre; 



2d 



dxdyK^^' (x, y) exp 



--{x-^ + y-O 



(C9) 



(CIO) 



is independent of e for K^^\x,y) strongly varying. In particular. 



i^^^)(e,0-(^/dxd,i.(f,f)exp 



--{x-^ + y-O 



i^(e,o- (cii) 
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Appendix D: Re^vriting the dissipator 

We shall derive the representation (43) of the dissipator. In what follows we shall drop 
the explicit notation of the momentum argument C, as well as the tildes in order to simplify 
the expressions. We note that we can express the matrix K with the help of the Pauli base 
(10) as 

K = Ko®l + K0a, (Dl) 



where K a = Kj ctj where i runs from 1 to 3 and we defined the components Kj G J^{ 
for i = 0, . . . , 3 according to Eq. (10). We express the remaining hermitian matrices Kj in 
a similar fashion in order to obtain 

Ki = kit + ki- a = fejol + hjCTj, (D2) 

and, hence, 

K = kool 1 + kojaj 1 + fcjol o"i + hjaj (Tj. (D3) 

It follows from the definition (10) that ky G M. Furthermore, from the indistinguishability 
of bath particles we note that Haa'piB' = i^a'ai3'i3 and, therefore, 

kij = kji i,j = 0,...,3. (D4) 

In a similar fashion we decompose the interaction potential V G t^(C)0e^(C) with respect 
to the second particle with the help of Eq. (10) as Eq. (15) where Vi G ,^(C) are given by 
Eqs. (15). With the help of these definitions we rewrite a typical sum which appears in Eq. 
(42) as 

1^1313' aa'VapVa' 13' = kijViVj, (D5) 

where the sum goes over all i,j = 0,...,3. We note that we can understand the above sum 
Eq. (D5) as the scalar product between a vector V = (Vb, Vi, V2, VJj)^ G ]R^0^(C) and the 
rotated vector /CV, where the matrix /C = {kij} G R^^^ is symmetric due to Eq. (D4), i.e. 
/C^ = /C. Since all elements of /C are real and since /C is symmetric, it follows that the /C 
may be diagonalized by an orthogonal matrix U G M^^^, where U^ = U^^ . Hence, denoting 
by (■, ■) the scalar product in M^, we have 

kijViVj = (V, /CV) = (V, W^T^WV) = (WV, TZUV) = (5, 7^5), (D6) 
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where we defined S = UV e M'^ ® ^(C). Denoting by Si the components of the vector 5, 
and by pi the diagonal elements of TZ where Si G ^^(C) and pi G M, we obtain 

l^p^'aaiyapVa'P' = PiSiSi- (D7) 

We remark that this convenient form of the sum Eq. (D5) is a resuh of the indistinguisha- 
bility of bath particles, which assures that K is a real, symmetric matrix. 

(2) 

Appendix E: The kernel of Q^ 

We shall briefly demonstrate that Q\ {A) = is equivalent to [5'i(^), A] = for all ^ and 
all i if detailed balance is required for A, where A, Si{^) G ^(C). Detailed balance means 
that each term contributing to Q^ , Eq. (46), vanishes individually. Hence, it suffices to 
investigate the operator 

Qf(A) = iy'de[[A(0,^],A(0], (El) 

where A(,^) G J^{C). We assume that A j^ 1 however Q)^ (A) = 0. Then, clearly 



tr 



Qf {A)a\ = I d^tr [AiOAAiOA - A(0 A(0^'] = 0- (E2) 

We rewrite the first term of this equation for all ^ as 

tr [A(OAA(OA] = A,,(OAfefc(0|a.fc(OI', (E3) 

where the Xa are the eigenvalues of A(^) and the aik{C,) are matrix elements of A represented 
in the eigenbasis of A(,^). In a similar fashion, we obtain for all ^ 

tr[A{OHOA']=XUOW^k{0\'- (E4) 

Since A(0 G ^2(C), 

tr [AiOAAiOA - A(OA(0^'] = " [Aii(0 - A22(0]' WMOl' < 0, (E5) 

and, therefore, 

Qf(A)<0, Qf\A) = 0^m),A] = Ve (E6) 

The statement [A(^), A] = ^ QJ^ '{A) = is trivial and, therefore, 

Qf\A) = ^ [5,(0, A] = Ve and Vz, (E7) 

if detailed balance is required for A. 
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Appendix F: Derivation of the Moyal product 



In this Appendix we briefly present the derivation of Eq. (74), i.e. the Moyal product, 
where the Moyal bracket has been specified in App. B. We start with Eq. (B3). We rewrite 
this equation in dimensionless variables, see Eq. (69), in order to obtain 



>V. 



Ap) {x,y) 



2d 



dx'dz' / d^d^a^^^ z' + -x',e + 7^ 



(^) .' 



1/7^ , 1 



(2vre; 
xiy(^)(^^'-la;',|-^e)exp 



'-{X - z') ■ i + '-X' ■ {^ - T]) 



^ I dx'dz' I dMa^^^ [z' + y, e + |e 



xiy(^) [z' -jx',i~ji]exp 



i{x- z')-i + tx' -{^-T]) (Fl) 



The expression (Fl) can now be expanded in a Taylor series in terms of e around e = 
under the assumption that o*-^-* and W^'^^ are both slowly varying functions with amplitudes 
of order one. Then, we obtain in zeroth order 



0, = a^'^W^'\ 



(F2) 



where Eqs. (C5) and (C6) have been used in order to eliminate of the integrals in Eq. (Fl) 
for e = 0. A similar calculation in first order demonstrates that 



O, = ^V.a(°) ■ V.l^W - ^V,aW ■ V.W'-'\ 



(F3) 



Higher order contributions can be obtained in a similar fashion but will not be discussed 
here. 

The same calculation can be carried out for the other term of the commutator and one 
obtains 



W, 



-V.H^(°)-V,a(°)]+0(£2). 



(0) 



(F4) 



Combining Eqs. (F2), (F3) and (F4) finally gives the desired result, Eq. (74), 

W, [C{A) (x,y)] = [a(°),iy(°)] + I [{a(°),iy(°)}^_^ - {W^^^U^^^} J + 0{e'), (F5) 
where {■, ■}x,r] denotes Poisson's bracket. 
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Appendix G: Derivation of Eqs. (80) 



We rescale Eq. (41) in order to obtain 



X 



Vpa (0 W^'^ ( X, r^ - ^, t ) V>,,, {O 



\%. (0 %'.' {O W^^^ U V - ^, t 



._]4/(£) j a;^7y 



^ + e' 



,tK3.(0Va'(n 



X exp 



:x ■ (e + e') 



(Gl) 



According to Eq. (C9), where we used that Vj^aiO is independent of e since V'^^^(r) is 
strongly varying in position space. Furthermore, from the definition of K, Eq. (25), we 
obtain that 

K^'\z, z') = ^C^'\z - z') + 6D(^, z'), (G2) 

where we already inserted assumption (78) and employed that due to Eq. (72c) D(^)(z, z') = 
bD{z, z'). Hence, with the help of Eq. (Cll) we obtain for the matrix elements of Eq. (G2) 

^(^)(e, o = rc^l.Aom + o + edl:l,,, (e, o ■ (G3) 

We shall now discuss these two contributions to the integral in Eq. (Gl) separately. We note 
that the Fourier transform of the strongly varying part of C*-^-', see Eq. (78), is independent 



of e, i.e. r, 



(e) 



Tq. Hence, if 7 = e the strongly varying part results in a collision integral 



of the form (80), for 7 = e^ or even weaker, this contribution vanishes and for 7 = 1 or 
stronger, the collision integral diverges. In a similar fashion, for the slowly varying part of 
C'-'^-* we have 



ri^HO 



dzTs{z) exp I — ^ ■ z 






(G4) 



hence, the resulting collision integral vanishes as e approaches zero. 

We shall now study the contribution to (Gl) arising from D. Applying the Fourier 
transform (Cll) to D [see Eq. (21)] gives the matrix elements 



^aa'BB'i^''^ ) 



^aa'^l3p' 



dznap{z) exp 



(27r£)2'i 



--z.{i + 



(G5) 
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Inserting this expression into Eq. (Gl) yields that this contribution tends to zero after the 

substitution C, ^ eC, and 

xi' — ;■ eC,'- The final form of Eq. (80) is obtained by performing the steps outlined in App. 

D. 

Appendix H: Derivation of Eqs. (84) 

In order to derive this result we rewrite Eq. (36) in rescaled variables as 

Vat3 (x + y - z\ w^^yy^,^, (x - y - Z' 



-\ye.p (x + ]^X' - ^ V^,p, (x + ]<t' - Z'\ W"^^' 
- V(^)V„;3 ix - ^X' - ^ V^,fi, [x - ]<t' - z^ 



X exp 



I 



x' ■ (jj — T]') 



(HI) 



e 

where W^'^' = W^^\x,'r]',t). We note that Va/^ir) is independent of e since V is slowly 
varying in position space. 

We now study the different scenarios arising from Eq. (82). For the slowly varying part 
of C*-'^'' we substitute x' — ?■ ex', draw the limit and integrate with respect to x' and rj' in 
order to obtain Eq. (84a) for 7 = e in the low density limit and for 7 = 1/e in the weak 
coupling limit. 

For the strongly varying part C-^-* we also substitute z ^ ez and z' — ?■ ez'. Hence, z 
and z' only appear in the matrix elements of Fq as e — ?■ 0. However, according to Eq. (26) 
this integral vanishes at all scales, thus, no contribution arises from the strongly varying 
part Fq in the semiclassical limit. Finally, in the low density limit the contribution of 
D(^)(2, z') = eJD^z, z') is easily seen to be of the form (84b). We remark this term vanishes 
in the weak-coupling limit. 
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